Spectral analysis of deformed random networks 
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We study spectral behavior of sparsely connected random networks under the random matrix 
framework. Sub-networks without any connection among them form a network having perfect com- 
munity structure. As connections among the sub-networks are introduced, the spacing distribution 
shows a transition from the Poisson statistics to the Gaussian orthogonal ensemble statistics of 
random matrix theory. The eigenvalue density distribution shows a transition to the Wigner's semi- 
circular behavior for a completely deformed network. The range for which spectral rigidity, measured 
by the Dyson-Mehta A3 statistics, follows the Gaussian orthogonal ensemble statistics depends upon 
the deformation of the network from the perfect community structure. The spacing distribution is 
particularly useful to track very slight deformations of the network from a perfect community struc- 
ture, whereas the density distribution and the A3 statistics remain identical to the undeformed 
network. On the other hand the A3 statistics is useful for the larger deformation strengths. Finally, 
we analyze the spectrum of a protein-protein interaction network for Helicobacter, and compare the 
spectral behavior with those of the model networks. 
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PACS numbers: 89.75.Hc,89.90.+n 



I. INTRODUCTION 

The network concept has been gaining recognition as a 
fundamental tool in understanding the dynamical behav- 
ior and the response of real systems from different fields 
such as biology, social systems, technological systems. 
Examples of biological systems include food-web, ner- 
vous system, cellular metabolism, protein-protein inter- 
action network, gene regulatory networks; social systems 
include scientific collaboration, citation, linguistic net- 
works, and technological systems include internet, power- 
grid PJ . Many of these networks have been shown to have 
universal structural properties, such as degree distribu- 
tion following a power law, small diameter, large cluster- 
ing coefficient, existence of communities Q], [2l H 0] ■ 

Different network models have been proposed and in- 
vestigated in detail to understand systems having an un- 
derlying network structure [ESSE]- These models 
concentrate to capture one or more structural proper- 
ties of the networks mentioned above [l|, ; @j- Apart 
from these direct measurements of structural properties, 
network spectra are also useful to understand various 
properties of the underlying system. Eigenvalues of the 
adjacency matrix of networks form what are called net- 
work spectra, and provide information about some basic 
topological properties of the underlying network |g, q}- 
Recently, considerable research has been done in the di- 
rection of network spectra [8, 9]. 

In the following, we mention known results on the 
spectra of real world and model networks. The spec- 
tra of networks have some correspondence with the spec- 
tra of random matrices. For instance, the distribution 
of eigenvalues of a matrix having finite mean number 
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p of nonzero Gaussian distributed random elements per 
row follows Wigner semicircular law in the limit p — » N, 
where N is the dimension of matrix [IOL llll ] . For very 
small p, which corresponds to the sparse random ma- 
trix, one gets the semicircular law but with peaks at 
different parts of the spectrum (maximum at the eigen- 
value zero) [l2(]- Recent investigations of the spectral 
behavior of networks, leading to matrices with entries 
zero and one, show that the random networks [131 ] fol- 
low Wigner semicircular law as well |14| with degener- 
acy at the eigenvalue zero. The small-world model net- 
works |2J show a very complex spectral density with many 
sharp peaks [Tj| , while the spectral density of the scale- 
free model networks (3[ exhibits a triangular distribution 
d, Q EM • The spectra of real world networks show 
remarkably different features than that of the model net- 
works [3, UM, [l6|, [13] , and based on this observation a net- 
work construction method was proposed which captures 
a peak at zero property shown by the spectra of many 
real world networks such as protein-protein interaction 
networks [l7| . Recently, spacing distributions of Erdos- 
Renyi networks have been studied under random matrix 
theory (RMT) framework [18|]. As connection probabil- 
ity decreases Ref. [18[ shows a transition to the Poisson 
statistics. Additionally, it shows the transition to the 
Poisson statistics upon the deletion of nodes in the real 
world networks 18]. Refs. [la LL9| have shown that the 
spacing distributions of various model networks, namely 
small- world and scalefree networks, follow the univer- 
sal behavior of RMT. In contrast to [la ], these works 
[la , [l3 ] have considered only connected networks. Fur- 
thermore, spectral rigidity such as the A3 statistics, de- 
fined in Eq. [3l provides a qualitative measure of the level 
of randomness in networks [201 ]. Recently localization 
of eigenvectors have also been used to analyze various 
structural and dynamical properties of real and model 
networks [211 ] . 

RMT, initially proposed to explain statistical proper- 



ties of nuclear spectra, has also provided successful pre- 
dictions for the spectral properties of different complex 
systems such as disordered systems, quantum chaotic sys- 
tems and large complex atoms among this. It has been 
followed by numerical and experimental verifications in 
the last few decades piJ, LLJJ- Quantum graphs, which 
model the systems of interest in quantum chemistry, solid 
state physics and transmission of waves, have also been 
studied under the RMT framework [13|. Recently, RMT 
has been shown to be useful in understanding the sta- 
tistical properties of empirical cross-correlation matrices 
appearing in the study of multivariate time series in sev- 
eral problems : price fluctuations in stock market [23J ]. 
Electro encephalogram data [2J|, variation of different 
atmospheric parameters [25J. 

In the present work we study spectral behavior of net- 
works having community structure under the framework 
of RMT. The study of community structure helps to 
elucidate the organization of networks, and eventually 
could be related to the functionality of groups of nodes 
d, |jj [26| . Regardless of the type of real world networks 
in terms of the degree and other structural properties 
1], it is possible to distinguish communities in the whole 
networks (4J. However, the question of definition of the 
community is problematic, and usually community is as- 
signed to the nodes which are connected densely among 
themselves, and are only sparsely connected with other 
nodes outside the community. We therefore model here 
community structure by sparsely connected Erdos-Rcnyi 
random networks. This simple approach considers more 
densely connected nodes as a definition of community, 
and does not pay attention to the detailed structure of 
the connections Q. Recent literature is largely filled up 
with methods to detect communities in networks based 
on structural measures [271 . |28| ] , whereas few works em- 
phasize on the spectral properties such as density dis- 
tribution and eigenvector analysis as well [29|]. The ob- 
jective of our work is not the detection of communities, 
rather we show the applicability of spectral methods un- 
der the RMT framework to analyze community struc- 
tures in networks. Instead of paying attention to the 
nodes forming communities, we look for the signatures 
of overlapping of communities in the spectra of the cor- 
responding adjacency matrix. We study various spectral 
behaviors, namely density distribution, nearest neighbor 
spacing distribution (NNSD) and spectral rigidity for de- 
formed random networks. We find that the NNSD de- 
tects even the small mixing of communities in the net- 
work, whereas spectral rigidity probed by the A3 statis- 
tics is suitable to analyze larger mixing, which is, in gen- 
eral, the case for real world networks. Communities are 
modeled by random or scale- free sub-networks, and inter- 
actions between communities are considered as random. 
For small interaction strength the NNSD of the network 
shows the transition from the Poisson to the Gaussian 
Orthogonal Ensemble (GOE) statistics. For large inter- 
actions, the A3 statistics shows systematic increase in 
the range for which it follows GOE statistics. Finally, 
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FIG. 1: (color online) Connection matrices corresponding to 
p = 0.01 and different values of q. (a) plots the connection 
matrix of the two sub-networks which do not have any con- 
nection between them, (b) corresponds to q/p — 0.1, and 
(c) depicts the case q/p — 0.5, when the connections between 
the the sub-networks are as large as 50% of the connections 
inside. 



as an application, we study the spectral properties of a 
protein-protein interaction network of Helicobacter un- 
der the RMT framework. 



II. DEFORMED NETWORKS 

For an unweighted network, the adjacency matrix is 
defined in the following way : Aij = 1, if i and j nodes 
are connected and zero otherwise. For undirected net- 
works, this matrix is symmetric and consequently has 
real eigenvalues. Random matrices corresponding to un- 
weighted random networks have entries and 1, where 
number of l's in a row follows a Gaussian distribution 
with mean p and variance p{l — p). This type of matrix 
is very well studied within the RMT framework [lCj, LL2J • 
We then turn our attention to the following structure: 
(1) Take to random networks with connection probabil- 
ity p\ the spectral behavior of the matrix corresponding 
to each of these sub-networks (blocks) separately follows 
GOE statistics. The matrix corresponding to the full 
network would be a to block diagonal matrix. (2) In- 
troduce random connections among these sub-networks 
with probability q. This configuration leads to m block 
matrix, with blocks having entries one with portability p, 
and off diagonal blocks having entries one with probabil- 
ity q. The above networks can be casted in the following 
form: 



A = A + A Q 



(1) 



Aq is a to blocks diagonal random matrix, where each 
block represents one community, and the off-diagonal 
block matrix A q denotes the interactions among the com- 
munities. Each block in A q is a random matrix, which 
for large N has mean q and deviation q(l — q). Since the 
nonzero values of q introduce deformation to the com- 
plete block diagonal form, we refer A being a deformed 
network. This terminology is motivated by the litera- 
ture on deformed random matrices [3fj| . Fig. [T] shows 
the connection matrices for m = 2 and various values 
of q. Fig. HJa) represents the two random sub-networks, 
each of size N = 500, with the connection probability 



inside a sub-network being p — 0.01 and between the 
sub- networks being q = 0. The ratio q/p, which can be 
considered as the relative strength of A q and Aq, mea- 
sures the deformation from the block-diagonal form of 
the matrix, or from the perfect structured network. The 
value q/p = 1, which corresponds to equal strength of 
inter and intra-community connections, yields complete 
random network. Fig. []Jb) plots the connection matrix 
for q/p = 0.1, which implies that inter-community con- 
nections are 10% of the intra-community connections. 
Fig. [IJc) shows the connection matrix for q = 0.005; 
for this value of q, the inter-community strength is 50% 
(q/p = 0.5) of the intra- community strength. Note that 
in numerical simulations we use the value of p equal to 
0.01, which leads to a sparse connected random network 
(N c ~ N) with the average degree < k >~ Nxp — 5, N c 
being the number of connections in the network. Larger 
value of p would lead to networks with the larger average 
degree. Real world networks are sparse [l| , and hence we 
chose such a small value of p. 



III. NUMERICAL SIMULATION RESULTS 

We denote the eigenvalues of the network by Xi, i = 
1, ... ,?7i x N, where N is the size of the sub-network, 
and m is the number of the sub-networks. Note that 
the size of each sub-network may be different, but for 
simplicity we consider here equal size. Fig. [3 plots the 
spectral density for m = 2 block matrices having qN 2 
non-zero off diagonal entries, corresponding to the two 
sub-networks connected with probability q. As discussed 
earlier q varies from q = 0, which corresponds to the 
two completely disconnected sub- networks ( A = Aq, 
Fig. [IJa) ), to q = p leading to a single random network. 
The cases for < q « p correspond to the configu- 
rations when the initial community structure is almost 
preserved. Increase in the value of q leads more entries 
of one in the matrix A q (Eq. [1} . Finally the q = p case 
destroys the community structure completely, and the 
network can be treated as one single random network. 
Fig. [2] presents the density distribution of eigenvalues for 
various values of q. The eigenvalues are scaled with re- 
spect to the spectra of the network for q/p = 1. With 
this scaling, the density distributions are not semicircu- 
lar for values of q < p. As the coupling between the two 
blocks increases (q > 0), the density distribution shows 
a transition to the semicircular form at q = p: 



pW = rW(^- A2 ). 



where Ao is the radius of the semicircular distribution for 
q = p calculated from the spectra of network as Aq = 



(A„ 



A mm )/2, X max and X min being the highest and 
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FIG. 2: (color online) Density distribution of the two random 
sub-networks connecting each other with probability (a) q — 
and q/p = 0; (b) q = 0.001 and hence q/p = 0.1; (c) q = 0.005, 
hence q/p = 0.5 and (d) q/p = 1 which corresponds to q — 
0.01 . Each block (random network) has size N = 500. The 
axes are scaled in such a way that the semicircle corresponding 
to q — p has unit radius (see text). All graphs are plotted for 
20 realizations of random sets of connections among the two 
sub-networks. 



A. Nearest neighbor spacing distribution 

In the following, we study spectral fluctuations of the 
networks for different values of q. In order to get univer- 
sal properties of the eigenvalue fluctuations, one has to 
remove the spurious effects due to variations of the spec- 
tral density and to work at the constant spectral density 
on the average. Thereby, it is customary in RMT to 
unfold the eigenvalues by a transformation A^ = N(Xi), 

where N(X) = J, p(A') dX' is the averaged integrated 
eigenvalue density [lOj. Unfolding is a transformation 
which produces the eigenvalues with a constant average 
level density. Since we do not have an analytical form 
for N, we numerically unfold the spectrum by polyno- 
mial curve fitting. 

Using the unfolded eigenvalues, we calculate the NNSD 
P(s), where sW = A^+i — A;, for different q values. Fig. [3] 
plots the spacing distribution for the two values of q, 
q = and q = 10~ 4 . For such small values of q, although 
the density distributions remain unchanged, the NNSD 
shows significant changes. Spacing distributions calcu- 
lated from the network spectra are fitted using Brody 
formula 13111 . 



j p /3 ( s )=As /3 exp(-as /3+1 ), 



(2) 



where A and a are determined by the parameter /3 as 
follows : 



the lowest eigenvalue. Now we turn our attention to the 
statistics of eigenvalue fluctuations. 



A = (1 + /j)aanda 
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FIG. 3: (color online) Nearest neighbor spacing distribu- 
tion for the two such values of q which gives two extreme 
statistics. Histograms correspond to the numerical values 
q = 10~ 5 {q/p = 0.001) and q > 10" 4 (<j/p > 0.01). The 
solid and dotted curves are, respectively, Poisson and GOE 
predictions of RMT. The figure is plotted for an average over 
20 realizations of the random set of connections between the 
networks. 



Eq.[5]is a semi-empirical formula characterized by the sin- 
gle parameter j3. As (5 goes from zero to one, the Brody 
formula smoothly changes from Poisson to GOE. As can 
be seen from Fig. [31 for q/p ~ 0.001(g ~ 10 -5 ), the value 
of the Brody parameter (3 ~ 0.2, which suggests that dis- 
tribution is very close to the Poisson [P(s) = exp(— s)] 
denoted by the dotted curve in the figure. As the value 
of q increases, j3 also increases, and it is of the order of 
1 for the value of q/p ~ 0.01 (which corresponds to the 
value of q as less as 10 ~ 4 ), and becomes insensitive for a 
further increase in q. For larger values of q, we analyze 
the spectra using the spectral rigidity test of RMT. 



B. Spectral rigidity via A3 statistics 

The spectral rigidity, measured by A3 statistics of 
RMT, gives information about the long-range correla- 
tions among the eigenvalues. The A3 statistics measures 
the least-square deviation of the spectral staircase func- 
tion representing the cumulative density N(X) from the 
best straight line fitting for a finite interval L of the spec- 
trum, i.e., 



A 3 (L;x) 



1 

— min 

L C1,C2 



x+L 



[N(X) - ciA - 02] d\ (3) 



where c\ and C2 are obtained from a least-square fit. Av- 
erage over several choices of x gives the spectral rigidity 
A 3 (L). For the uncorrelated eigenvalues, A 3 (i) = L/15, 
reflecting strong fluctuations around the spectral density 



FIG. 4: (color online) Long-range correlations among eigen- 
values. Different open symbols are the numerical values of 
A3 for various q values, and the solid curve (merged with the 
open circles corresponding to q/p = 1) is the GOE prediction 
(Eq. [3} . Since for q — p the A3 statistics of network follows 
the GOE prediction completely, the solid line showing GOE 
statistics merges with the circles showing numerical values for 
this q. The figure is plotted for an average over 20 realizations 
of the networks. A3 follows the universal RMT prediction up 
to certain L values. The range of L for which A3 follows GOE 
statistics increases with the ratio q/p. 



p(A). For the GOE case, As(L) statistics is given by 

1 



A 3 (L)~ — InL. 

■K 1 



(4) 



Fig. Q] plots the A 3 statistics for five different values of 
q. Various open symbols are the numerical values of A3 
for various q values, and the solid line (merged with the 
q/p = 1 case) is the A 3 (L) statistics for the GOE case 
(Eq. 15). As seen from Fig. HI the A 3 (L) statistics follows 
RMT predictions of GOE (Eq. g]) up to a certain L. It 
has a linear behavior in semi-logarithmic scale with the 
slope of ~ I/71" 2 . The value of L for which it follows GOE 
statistics depends upon q. For small values of q such as 
q/p = 0.01 and q/p = 0.05, A 3 follows RMT prediction 
till very small range of L ~ 5 and L ~ 20, respectively. 
As q increases, the value of L for which A 3 follows the 
GOE statistics also increases. For q/p = 0.1, it agrees 
with the RMT predictions of GOE behavior for L ~ 75, 
and after this value, deviation from the RMT prediction 
is seen. This deviation corresponds to the existence of 
community structure in the network. As the value of q 
increases, the communities have more and more random 
connections between them. For q = p the community 
structure is destroyed fully, and the network is a com- 
plete random network. This fact is reflected in the A3 
statistics corresponding to q/p = 1. At this value of q, it 
follows RMT prediction up to a very long-range L ~ 150. 
After this value of L, for the network of size N x m = 1000 
we do not have a meaningful calculation of the A3 statis- 
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FIG. 5: (color online) Spectral behavior of deformed scale- 
free networks. Sub-networks are scalefree networks of size 
AT = 500 and average degree < k >= 5. (a) plots the eigen- 
value density distribution of deformed scalefree network for 
various values of q/p. (b) and (c) plot the NNSD and A3 
statistics respectively. Dashed and solid line in (b) corre- 
sponds to Poisson and GOE statistics respectively. All graphs 
are plotted for 20 sets of random realization of interaction net- 
works, (d) plots the density distribution (inset), and spacing 
distribution for a protein-protein interaction network in Heli- 
cobacter as histogram and GOE statistics as solid line. This 
network has size AT = 712 and average degree < k >~ 5. 
Open circles in (c) is the A3 statistics for Helicobacter. 



tics [32|. For q = 0.005 (q/p = 0.5) (see Fig. [TJ, where 
the strength of inter-community is as large as ~ 50% of 
the intra-community connections strength, the A3 statis- 
tics correctly reflects the deviation from complete ran- 
dom matrices, suggesting the existence of communities 
in the network. 

Note that we present results for each sub-network hav- 
ing equal size. For sub-networks having different sizes, 
all the figures remain the same. The crucial quantity, 
which affects correlations of eigenvalues, is the variance 
of each block or the ratio q/p. For blocks with different 
sizes, but with the same q/p, similar results are obtained, 
except for the exact value of L in Fig. 0] for which the A3 
statistics follows GOE distribution, which scales with the 
network size [20j . 



IV. DEFORMED SCALEFREE NETWORKS 

In the following we consider scalefree networks as the 
sub- networks, and study the spectral behavior for vari- 
ous values of q. Again q measures the strength of the off- 
diagonal block matrix defining the interaction between 
the sub-networks. Matrix Aq in Eq. [TJ corresponding to 
the scalefree sub-networks, consists of two block diago- 
nal matrices, with entries of one in each block follow- 



ing a power law characteristic of the sub-network. We 
use Barabasi- Albert algorithm 0] to generate the scale- 
free sub- networks. In scalefree network the probability 
P(k), that a node has degree k, decays as a power law 
P(k) ~ k ~ 7 , where 7 is a constant and for the type 
of probability law used in the simulations 7 = 3. Other 
forms for the probability law are also possible which gives 
different exponent 33]. However, the results reported 
here are independent of the value of 7 [3J| . Size and av- 
erage degree of the sub-networks remain the same as for 
the random sub-networks, i.e. N = 500 and < k >= 5. 
The average degree (< k >) of a network can be cal- 
culated as < k >— 2 x N c /N, where N c is the number 
of connections and A?" is the size of the network. With 
the increase in the value of q, deformation from the net- 
work having scalefree community structure also increases. 
Fig. [5] plots various spectral behavior of deformed net- 
works made of the scalefree sub-networks. Figf5]Ja) plots 
the density distribution for the various values of q. For 
small values of q, the density is very different from that of 
the deformed random networks (Fig. [2]). It has a triangu- 
lar shape with a peak at zero. This is a well-known shape 
for sparse scalefree networks [9|, [l5|, LL6|- For q/p < 0.01, 
when the scalefree structure of the sub-networks dom- 
inates over the random interaction between them, the 
eigenvalue density distribution does not show any no- 
ticeable change. But the NNSD in Fig. [5{b) suggests a 
possible structure in the network. As shown in Fig.[5jb), 
for q/p = 0.001 (q = 10~ 5 ) the NNSD is close to Poisson 
statistics with a value of the Brody parameter j3 ~ 0.21. 
As q increases, value of the Brody parameter increases 
as well, becoming one for q ~ 10~ 4 . After this value of 
q, the NNSD does not provide any further insight, and 
we probe for long-range correlations among eigenvalues. 
Fig-EI c ) plots the A 3 statistics for various values of q. It 
shows similar behavior as for the deformed random net- 
works (see Fig. @|. For q ~ 0.01, when the network has 
distinguishable community structure, the value of L for 
which A3 follows the GOE statistics (JH) is as small as 25. 
As q is increased, L also increases, becoming ~ 150 for 
q/p~ 1. 

Fig. [5[d) shows the density distribution (inset) and 
the spacing distribution of the protein-protein interaction 
network of Helicobacter [35| . The largest connected com- 
ponent of the network has dimension N = 708 and num- 
ber of connections N c — 2789. The average degree of this 
scalefree network is < k >~ 4. The density distribution 
has triangular form with a peak at zero. This behavior 
of the density distribution suggests scalefree properties 
of the network [9|, [la, [lg] , but does not provide informa- 
tion of randomness or structure in the network. To get 
further insight, we calculate the NNSD and the spectral 
rigidity of the network. For this, first we unfold the eigen- 
values using the procedure explained earlier. The NNSD 
of the network follows GOE statistics with the value of 
(3 ~ 0.98, suggesting enough random connections in the 
network. Further test of long-range correlations among 
eigenvalues shows that the A3 statistics follows the GOE 



prediction (Eq. [4]) up to L ~ 20 in Fig. 02c), and after 
this value deviation from the universal behavior is seen. 
It suggests that, though the network has enough random 
connections which give rise to short-range correlations 
among eigenvalues, it has strong community structure 
causing deviation of the A3 statistics from the random 
matrix behavior after a certain range. 

In the present paper we consider only the random in- 
teractions between communities. For other kind of in- 
teractions, for instance interactions among the scalefree 
sub-networks as considered in [5j which leads to a hierar- 
chical scalefree network, the density distribution would 
show an entirely different behavior from the semi-circular 
distribution, p = q case would lead to a scalefree topol- 
ogy which has a triangular density distribution with peak 
at zero. However, spectral fluctuations would show qual- 
itative similar behavior. For small coupling interactions 
among the sub-networks, the NNSD results would be 
same as presented here, showing a transition from Pois- 
son to GOE statistics [20|], whereas for large coupling 
interactions the exact range for which the A3 statistics 
follows GOE would be different from those of the ran- 
dom interactions. Further detailed results of this model 
as well as real world networks having more complicated 
structures analyzed under the deformed random matrix 
framework would be discussed elsewhere 1 36] . 



V. CONCLUSIONS AND DISCUSSIONS 

The eigenvalue density distribution of networks hav- 
ing two sub-networks tend towards the semi-circular dis- 
tribution as the random connections between the sub- 
networks are increased. For very small values of q < 
1CP 4 , corresponding to the very small deformation from 
the community structure, the density distribution does 
not present any noticeable changes, but the NNSD, which 
reflects short-range correlations among eigenvalues, show 
important features. For two random sub- networks, which 
are almost uncoupled (i.e. q ~ 0), the NNSD is very 
close to the Poisson statistics, and as q increases, it has 
a smooth transition to the GOE statistics. Note that 
this Poisson to GOE transition is found for many differ- 
ent systems, for example spectra of insulator-metal tran- 
sition, order-chaos transition follow this Poisson-GOE 
transition [ll|. Sade et. al [37J have studied transi- 
tion to the GOE statistics as a function of site disorder 
for the spectra of small-world and scale-free networks. 
Here, by keeping the network structure fixed, disorder at 
nodes is increased and depending upon the network av- 
erage degree transition to GOE statistics is seen. The 
main difference between 37} and the study presented 
in this paper is the following: we track changes in the 
spectra with structural changes in the network architec- 
ture. As random connections among the sub-networks 
are increased, first there is transition for the NNSD to 
the GOE statistics, and this transition occurs for very 
small value of random connections among networks. This 



is the crucial and remarkably different result observed 
here, which suggests that very small random interaction 
between communities is enough to introduce short-range 
correlations among them, spreading the randomness in 
the whole network. Second, further increase in coupling 
among the sub-networks is reflected by long-range corre- 
lations among eigenvalues. For this increase in the value 
of q, the NNSD does not give additional insight to the 
deformation of the network, as it remains same with the 
/? ~ 1, so we turned our attention to the A 3 statistics. 

The A3 statistics, which measures long-range correla- 
tions among the eigenvalues, detects deformation from 
a network having two coupled sub-networks, to a single 
random network. More deformation of the network from 
community structure, leads to a larger range of L for 
which A3 follows the GOE statistics. Note that, for the 
case of sub-networks being completely random, the spac- 
ing and the A3 statistics of each of them follows RMT 
prediction. Therefore, any deviation from GOE statistics 
is due to the community structure these two sub-networks 
form when considered as a single network. 

It is interesting to note that our results resemble the 
behavior of deformed random matrix ensembles (DGOE) 
introduced to study the effect of isospin symmetry break- 
ing in nuclei [301 ] . The qualitative behavior of the spec- 
tral density and the A 3 statistics of networks presented 
here is similar to that of deformed matrices studied in 
1381 KM |40| . The analytical form of the density derived 
in [391 ] depends on a parameter a measuring the rela- 
tive strength of the off-diagonal random matrices to the 
block diagonal random matrices. In similar lines, for de- 
formed networks, we can compare q/p, relative strength 
of off-diagonal and diagonal networks, with a. The re- 
sults presented here suggest that further investigations of 
complex networks following similar lines as in deformed 
random matrices [39j would be useful to have detailed 
information of communities in the networks [4l|] . 

To conclude, we have studied the spectral behavior of 
networks having community structure, and shown that 
the NNSD and A3 statistics capture features related to 
the structure in the network. We investigate the spectral 
properties of a real world network as well, and compare 
the results with those of the model networks. On the one 
hand, results presented in this paper advances the studies 
of the spectral properties of network with the community 
structure under the universal RMT framework; on the 
other hand, variations in the correlations among eigen- 
values shed light on the coupling among communities. 
For the simulations, the community structure in network 
is modeled by the very simple random or scalefree sub- 
networks, and the interactions among these sub-networks 
are considered random, whereas real world networks have 
richer structure [5j. However, the results presented here 
provide a platform to investigate the community struc- 
ture of networks using a well developed theory of ran- 
dom matrices; the further investigations in this direction 
would deal with real world networks with richer and more 
complicated structure under the deformed random ma- 



trix framework [39, [4l|| . 



prospective of the results. 
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